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J> ' One possible way in which phase transitions in the early universe may have ocurred is 

o\ : 

■ via nucleation of bubbles of the new phase (true vacuum) in the old phase (false vacuum). 

O ! The technique most widely used to compute the probability of bubble nucleation is based on 

^ . 

' instanton methods in the context of the semiclassical approximation. At zero temperature in 

\ 3 + 1 dimensions the nucleation rate is dominated by the 0(4) symmetric instanton, a sphere of 

radius R, while at temperatures T » R~^, the decay is dominated by a "cylindrical" (static) 

D , instanton wtih 0(3) invariance. There has been discussion in the literature as to whether the 

' transition between these two regimens would be first order (discontinuity in the first derivative 

• 1—1 , 

^ . of the nucleation rate at the transition temperature Tc), or second order (continuity of the 

' first derivative, but discontinuity of the second derivative at Tc). In this paper we obtain 

the finite temperature solutions corresponding to the quantum and the thermal regimes, and 
compute their action as a function of the temperature for different values of the wall thick- 
ness in a if^ potential. Our results indicate that only for the cases of very large wall thickness 
a second order transition takes place, while for all the other cases a first order transition occurs. 
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1 Introduction 



The problem of the decay of a metastable state via quantum tunnehng has important apphcations 
in many branches of physics, from condesed matter to particle physics or cosmology. In many cases 
however, the difficulties associated with a full quantum treatment, and the fact that tunneling is 
a non-perturbative effect, make it impossible for us to be able to make quantitaive predictions 
about the behaviour of the system. When the nature of the tunnehng system allows the use of 
the semiclassical approximation though, instanton techniques, used to describe quantum tunnel- 
ing in that regime, provide us with a formalism capable of producing accurate values for decay 
rates. In condesed matter this formalism has been applied to an array of problems ranging from 
the tunneling of an Abrikosov vortex line out of a pinning potential in a superconductor [1-3], to 
the change in the magnetization direction of thin films under field inversion [4] . In the context of 
cosmology, perhaps the most compelling example is that of the study of first order phase transi- 
tions in the early universe [5] , although other problems such as the decay of metastable topological 
defects formed during these transitions [6], have also shown to be amenable to this sort of treatment. 

The basic tenet of the formalism is that the decay rate per unit time (and where applicable per 
unit volume) is given by an expression of the form 

r = Ae-^^ (1) 

where Se is the Euclidean action of the instanton: the classical solution to the equations of motion 
under the potential barrier in euclidean time with periodic boundary conditions. The instanton 
has turning points at the configurations at which the system enters and exits the barrier, and its 
analytic continuation to Lorentzian time at the exit point describes the appearance of the system at 
the other side of the barrier and its subsequent evolution. When the solution for the given boundary 
conditions is not unique, one should take into account only the one with the lowest Se, which will 
dominate (1) (unless of course all or some of the solutions happen to have comparable actions, in 
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which case one should take all relevant contributions to F into account). The determinantal pref- 
actor A in Eq.(l), meanwhile, results from a gaussian functional integration over small fluctuations 
around the instanton. The formalism was first developed and applied to false vacuum decay at zero 
temperature in refs [7-9]. The dominant solution in this case is a maximally symmetric instanton 
invariant under 0(4) (for three spatial dimensions) [9]. 

In the context of quantum mechanics (i.e. + 1 dimensions) the formalism was extended to include 
finite temperature effects by Affleck [10]. Under certain assumptions for the shape of the barrier he 
found that above a critical temperature determined by the curvature of the potential at the top of 
the barrier, the dominant solution was a "static" (euclidean time independent) instanton sitting at 
the bottom of the inverted potential well, the decay rate produced by this solution recovering the 
expression for purely classical thermal activation over the barrier that we would expect to dominate 
the procccs at high temperatures. The transition between the thermal and quantum regimens is 
dominated by solutions with a finite period in the Euclidean time that smoothly interpolate between 
the zero temperature instanton and the "static" high temperature solution. The transition is thus 
a second order one (i.e., the decay rate F and its first derivative with respect to the temperature 
are continous at the critical temperature of the transition, Tc, but not so the second derivative). 

Chudnovsky [11] however showed that this situation is not generic, and that the type of transi- 
tion will in general be determined by dTp/ dE^udi the derivative of the euclidean period of the motion 
Tp with respect to its euclidean energy. For the case in which the period vs. energy curve possesses 
a minimum at values of {Eeuci,Tp) different from those associated with small oscillations around 
the minimum of the potential, it can be shown that, at the transition temperature, the dominant 
solution changes discontinously to the high temperature instanton. Consequently the curves for the 
actions of both regimens cross each other at that point, and the derivative of Smin with respect to 
the temperature (and thus that of F) changes discontinuosly as we go from the thermal activation 
regime to the quantum regime. This is what Chudnovsky called a first order transition. In this 
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case we should be able to observe quite a sharp change in the behaviour of F as we go through the 
critical region. 

In the context of field theory very little work has been done towards examining the structure 
of the intermediate regime instantons as well as that of the transition between the two extreme 
cases of high and zero temperature. Studying the bubble nucleation rate for phase transitions in 
the early universe, Linde [5] suggested that periodic instantons would smoothly interpolate between 
these two regimens, much in the same way that Affleck had shown could be the case for quantum 
mechanics shortly before. In this picture then, one would have Coleman's 0(4) invariant spherical 
bounce of radius R dominating at zero temperature. As one goes into finite temperature one would 
first encounter a periodic array of, rather weakly interacting, widely separated spheres, but after we 
have increased the temperature up to values T ~ {2R)~^ and beyond we no longer can fit the spheres 
into the corresponding periods, and they would start merging into one another producing what has 
become known as a "wiggly cylinder" solution. As we keep increasing the temperature further, the 
wiggles eventually smoothly straighten out, and the solution goes into the 0(3) invariant "cylinder" 
that dominates the thermal activation regime. 

In correspondence with the work of Chudnovsky in quantum mechanics however, Garriga [12], 
in a recently published paper, has shown that within the thin wall approximation the transition is 
in fact not second order, but rather, first order. The wiggly cylinder solutions do indeed exist within 
some range of the temperature, but when they exist they always have a higher action than any of 
the other solutions for the same temperature, and therefore they never dominate the exponential in 
(1). The instanton that saturates (1) goes then abruptly from being a periodic series of separated 
bounces to being the cylindrical solution, and the rate of change of the action with the temperature 
changes discontinuosly. The transition thus produced is in consequence first order. It may be worth 
noting at this point that the terminology used to describe this transitions (i.e., first and second 
order) refers only to the type of transition between the two different instanton regimes, and not to 
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the phase transition of the physical system itself. 

In this paper we aim to explore the nature of the transition between the thermal hopping and 
the zero temperature regime when we go beyond the thin wall approximation, in the context of 
the bubble nucleation problem in c/?^ models just mentioned. In order to do this, we will develop 
an algorithm to numerically compute the relevant finite temperature solutions for a given value 
of / (the parameter that governs the width of the wall as compared to the radius of the bubble), 
and then compute the corresponding curve for the action as a function of temperature. Finally we 
repeat this procedure for different values of /, covering the most significant physical cases. Our 
results indicate that in the case of extremely thick wall (when the tunneling occurs over a very 
shallow barrier) we do get a second order transition, but that in all other cases the transition is in 
fact first order, in particular, close to the thin wall regime our results agree rather well with those 
of Garriga. In spite of the fact that in any realistic model of the phase transition the potential will 
be temperature dependent, whereas we use temperature independent potentials, we do not expect 
this to significantly change our results, as will be seen at the conclusion. The paper is organized 
as follows: in section II we give a quantitative account of the finite temperature transition in the 
spirit of [11], in section III we present the results that Garriga obtained for bubble nucleation in 
the thin wall limit, in section IV we introduce our numeric algorithm and in section V we present 
our results. Finally, conclusions are in section VI 

2 Finite Temperature Instantons in Quantum Mechanics 

Acording to [7] the decay rate in the semiclassical limit {S » U) is, 

V(xexp{-Smin/Ti) (2) 
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where SminiT) is evaluated along the g(r) trajectory with period Tp = {3 = h/ksT that minimizes 
the Euclidean time action 



with q = dq/dr and where the integral is performed over a whole period (3 (sec [7-9], [11] for 
details; from now on we will use /? to designate both the euclidean period and the usual inverse 
temperature). Such trajectories will of course be solutions of the classical equations of motion in 
Euclidean time 

dV 

Mi^- (4) 

consistent with the boundary conditions, that is, periodicity with period p. (i.e., qi = g(— /9/2) = 
q{j3/2) and §2 — ?(0) will be the turning points of the motion). By noting that (4) is the equation 
of motion of a particle in a potential -V in Lorentzian time, we can draw two consequences: first, 
that the (positive) energy 



is conserved. And, second, that by inverting the barrier we can recognize the types of solutions 
that we will have for the motion of the particle in the well (see Fig. 1). It is then immediately 
obvious that there are two solutions with period /3, one that actually is a time independent solution 
sitting at the bottom of the well, q = qo, and another one that corresponds to periodic motion in 
the inverted well between qi{E) and q2{E). For the first one we will have from (3) an action 




(3) 




(5) 



5-0 = V{qo)l3 



(6) 



whereas for the second one, use of (5) into (3) yields 




(7) 



where E and (5 are related through 




(8) 
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Then, from (6), (7) and (8) one obtains 

dSo _ dSr 
d^S ~ ° ' 'dp 

and 

d'So _ fSr 
dp^ ' 

The behaviour of these derivatives at the temperature at which St = Sq will then determine the 
type of transition between the two regimens: if at Tc we have S'q^S^ the transition will be first 
order, but if the equality holds, we will have a second order transition for which S'^^S'^ as we are 
about to sec. From (10) however we see that this will in turn be determined by the behaviour of 
dj3 1 dE (which can be interpreted as the specific heat of the system) . 

As a first example we can look at potentials of the form — + (f" or —(^ + [10]. In this 
case /5 monotonically decreases with E (Fig. 2), until it eventually reaches /?o, the period of small 
oscillations at the bottom of the potential. As seen from Fig. 2, St—>-So and S^^Eq = S'q as P^Po, 
or what is the same, as T^Tq, Tq{— h/kBPo) being the transition temperature in this case. The 
second derivative S'^ will however change discontinually as we go across Tq, and will, as a matter 
of fact be undefined at that temperature (of course a more careful quantum analysis is needed at 
this point). This is what Chudnovsky called a second order transition. 

The other case of interest to us is that for which P{E) develops a minimum, at say {Ei,(3i), 
before reaching /?o at £^0 (see Fig. 3). This will be the case for potentials that change very slowly 
both at the top and the bottom (although still being parabohc) but that are very steep in the 
middle. In this case we see that, as E decreases from £"0 to Ei, we have = E < Eq = Sq and 

— 1/ P' > 0. The corresponding curve for St{P) will then be a concave curve that stays always 
above the thermal activation curve. At Ei however P' — and is undefined, and for E < Ei 
we have to have a convex curve that quickly goes into a straight horizontal line for high vlaues of 
P (Fig. 3b). St then has to cross over the Sq curve at some period Pc- The derivatives of the 
two curves will be different at that point and, therefore, the first derivative of Smin will suffer a 
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= £ > , 



dp/dE 



(9) 



(10) 



discontinoiis change as we go across the corresponding critical temperature Tc = h/ksf^c (note that 
Tq < Tc < Ti), producing thus a first order transition. 

3 Finite Temperature Instantons in Field Theory 

Studying first order phase transitions in the early universe, Linde [5] suggested an extension of the 
previously discussed second order transition scenario (published by Affleck shortly before) to the 
field theory case. Due to spherical symmetry in the spatial coordinates, the equation of motion for 
the field is now 

5V 2dip dV 

— — H — H — = — , (11) 

dr^ dr^ r Or dip ^ 

where r is the spherical radial distance. We will choose the potential to be of the form: 

V{^) = ^{ip'-,x'f-Fip . (12) 

As for the boundary conditions, in the r direction the requirement of finite action leads to ip^ipj^ 
as r^oo, while in the r direction periodicity alone no longer suffices, and we have to impose that 
the configurations atT = ±/3/2, have zero derivative along the r direction to ensure that they 
are indeed turning points. We thus have 

as r— >-oo , d(p/dT = at r = ±/3/2 , , (13) 

where j3 is the period of the solution. The idea is then that the relevant solutions should follow 
the same lines as above: at low temperatures the solution would be a periodic array of widely 
separated 0(4) spherical bounces (with radius R). As we go into temperatures T>{2R)~^ however 
we no longer can fit a bounce in each period and they start merging into one another, forming an 
oscillating "wiggly" cylinder analogous to the oscillating solutions of the quantum mechanical case. 
These wiggles would then subsequently smooth out as we increase the temperature and eventually 
disappear at temperatures T » R~^, providing a smooth transition into the "cylindrical" (r 
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independent) instanton corresponding to the thermal activation regime. Fig.4 shows this sequence. 
As we increase the temperature, the transition would be then a second order one. 

As noted in the introduction however, Garriga [12] has shown that at least in the thin wall 
approximation this is not the case. In this approximation, the parameter F in the potential above 
(and correspondingly the energy difference between the two different phases, e = 2F) is very small, 
and as a consequence the width of the wall separating the bubble of true vacuum from the false 
vacuum is very small compared to its radius R (i.e., one needs a larger bubble to compensate for 
the wall energy). In Lorentzian time then, the dynamics of the bubble can adequately be described 
by an effective action (in 3 spatial dimensions) 



(14) 



where the first term is the Nambu action proportional to the area of the world sheet of the wall (cr 
being the wall tension, a set of coordinates on the world sheet and 7 the determinant of the world 
sheet metric) , while the second is the volume of the bubble times the energy difference between the 
two phases integrated over time. With a fiat background geometry, adopting spherical coordinates 
and a spherical ansatz for the world sheet means that we can write S as 



S = -aSo 



(15) 



where f — dr/dt and S2 is the surface of the unit sphere. Since the Lagrangian does not depend on 
time, the energy 



E — Prf — L — (TS2 



'I _ ^2)1/2 2,a 
must be conserved. With a little algebra this can be written as 



(16) 



+ 1/(r, E) = Q , 



(17) 



where 



V 



E 
(J So 



- 1 



(18) 



The instantons relevant to the problem are then found by performing a Wick rotation in (17), 
whence we obtain 

-V(r,E) ^0 , (19) 

where now r = dr/dr. There are two types of solutions to (19) depending on the value of E: for 
E — and /3— >oo, the solution is just Coleman's bounce with radius R — 3a/e, and action 

ST^aS-iR^/A^a7r^Ry2 . (20) 

Keeping E — but introducing a finite value for /3 produces then a periodic array of bounces in the 
T direction. These solutions exist for as long as T<T* = {2R)~^ and will continue to have the same 
action St, since outside R the field decays exponentially. For E > though, the inverted potential 
-V takes the form of a well and the motion of r in it acquires two turning points, rmin and Vmax- 
The solution is then wiggly cylinder shaped. If we keep increasing E, —V becomes shallower and 
shallower, and the two turning points start approaching each other until they eventually merge at 
values EqjTq. From the condition V{rQ,EQ) = V'{rQ,EQ) = {V = dV/dr), we get values for the 
radius of the cylinder 

ro = ^i? , (21) 

and its corresponding energy 

£i = l5l2!<.= ^.i;^ . (22) 

Since it is static, the solution r — rg will exist at all temperatures, and its Euclidean action 

So = Eo/3 (23) 

will be proportional to the period (= temperature^^). This will therefore be the dominant solution 
at high temperatures. 

It suffices now to compare St and 5*0 to see that there will be a temperature for which St — Sq, 
namely 

32 

r. = ^«-' . (24) 
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At this temperature then, the curves for 5*0 and St will cross each other, and consequently the 
slope of the curve that follows the minimum action Smin will change discontinously. The transition 
between the two regimes is thus first order. It can be shown (see Ref.l2) that the wiggly cylinder 
solutions that exist for < ii^ < Eq, do, as a matter of fact, always have greater actions than either 
the periodic array of bounces or the static cylinder, and therefore never dominate Smin- 

The question that we aim to answer is then the following: Is this situation general for bubble 
nucleation, or will we find a second order transition as suggested by Linde when we go into the 
thick wall case?. 



4 Algorithm 



We start with the Euclidean action 



S, = A.jdrjdrr^ 



where if = dip/dr, (p' = dip /dr. We then rescale the Lagrangian 



(25) 



(26) 



to get 



(27) 



The only adjustable parameter in the Lagrangian is now /, so by covering its whole range we should 
be covering all relevant cases. The equations of motion are now 

5V <9V '2 dip dV 



dr'^ '''' r dr dip 
and the boundary conditions are the usual ones 



2{<p'-<p)-f , 



(28) 



</7^</7+ as r— >oo , d(p/dT — at r = ±/3/2, 



(29) 
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To solve this equation we have used a multigrid algorithm. Multigrid methods were first introduced 
in 1970s by Brandt [13] (see also [14] for a brief introduction and further bibliography on the 
subject), and the basic underlying idea is the following: let's suppose that we have to solve a linear 
equation of the form 

Lh{qh)^0 , (30) 

where L is some discrete form of the linear differential operator, q is the discrete solution, and the 
subscript h denotes a discretization of the independent variables (i.e., r and r in our case) in a mesh 
with meshsize h. In order to do this, iterative methods usually start with an approximate solution 
Qh such that 

LhiQh) = dh , (31) 
where dh is called the defect. We then look for approximate solutions, Vh, of 

Lhivh)^-dh , (32) 

to form 

qr-Qh + Vh . (33) 
quew then be closer to the true discrete solution than what Qh was. Under some general 
conditions for the form of the matrix of coefficients of the discretized algebraic equations, iteration 
of this procedure is a convergent process, converging upon qh- The heart of the matter is then what 
approximation to use to solve (32) and find Vh- 

As an example we have relaxation methods, where what one does is to replace Lh in (32) by a 
simpler operator (such as just its diagonal part for instance), and then solve the set of equations. 
The drawback here however is that, as solving the set of discretized equations is basically a point 
by point procedure, convergence for the long wavelenght modes of the solution will be rather slow 
usually. 

Multigrid methods operate by coarsifying rather than simplifying the problem. Instead of solving 
the problem in a mesh of size h with a simpler operator, we project the problem into a coarser grid 
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of size H, solve there for vh and then interpolate back to find the corrections Vh in the finer grid. 
Appart from having to solve far less equations than in the finer grid, transfering to a coarser grid 
has the added advantage of dramatically improving the convergence rate of the long wavelenght 
modes. However, if applied as such, this process will not only not converge, but slightly diverge for 
short wavelength modes. What one does then is to apply one or more realaxation sweeps before 
and after transfering to the coarse grid, to smooth out the short wavelength componenets of the 
error that relaxation handles so well. 

Several developments on this basic scheme are now possible. In first place, it is obvious that 
nothing prevents us from adding more grids to the algorythm, going each time to coarser and coarser 
meshes. In principle nothing prevents us in fact from going to grids so coarse that we only have to 
solve a handful of equations. Starting from such a grid and computing in it not just the correction, 
but the full solution qh, and then going up the ladder and so on, eliminates the need for having an 
initial guessed solution fed into the algorythm. This is known as the full multigrid method. 

Another possible expansion leads to an algorythm that can solve non-linear problems. In this 
case, if we have to solve 

Lh{qh) = fh , (34) 

where fh is some right-hand side term, we need a smooth correction Vh to the current solution 

so that 

Lh{qh + Vh)^fh ■ (35) 

To do this note that 

Lh{qh + Vh) - Lh{qh) ^ fh - Lh ^ -dh ■ (36) 
Transfering this equation to the coarser grid after the pre-relaxation sweeps leads to 

LhM - LniJlUh) = -nUh , (37) 
where TV^j is the restriction operator from gridsize h to gridsizc H . That is, in the coarse grid we 
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solve for 



(38) 



If the approximate solution is then qh, then the coarse grid correction is 



(39) 



and the corrected fine grid solution will then be given by 



Qh — <lh+ -Lh vh 



(40) 



where X// is now the interpolation operator between gridsize H and gridsize h. Note that in general 



For our problem we used full weighting injection for the projection operator (see Ref.l4) and second 
order polinomial interpolation (fourth order for the thin wall case). We were precluded from using 
a full multigrid approach because we encountered a bifurcation problem. This problem can appear 
in cases where the non-linear equations have more than one solution compatible with the boundary 
conditions, and in practice it means that for certain values of the mesh size h the coefficients of the 
discretized equations blow up. One is therefore confined to using meshsizes smaller than that for 
which the problem appears, and does not have acces to a mesh coarse enough to be able to solve 
the equations exactly. To solve this problem we limited our algorythm to use only three meshes. 
As for the needed initial guess solution, what we did was to obtain the zero temperature profile of 
the bounce via a simple shooting routine (remember that at zero temperature the dominant solu- 
tion is 0(4) symmetryc, and therefore the problem can be reduced to a 1-dimensional differential 
equation). We then put it in a two-dimensional form and fed it to the multigrid algorythm. After 
computing the form of the zero temperature solution in this way we then proceeded to compute 
finite temperature solutions by continuation in the temperature (i.e., using the previous solution as 
the guess for the next one with slightly higher temperature). 



in^l, and so qJl^"^ ^ IhQn- 
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Since Zebra relaxation turned out to be unstable for our problem, we used line Newton- Jacobi in 
the T direction for the relaxation smoother (as increasing the temperature meant in practice reducing 
the mesh size in that direction) . Special care had to be taken to handle the negative modes of the 
instantons. One of the general conditions that relaxation methods need to converge is the absence 
of eigenvalues with opposite signs in the diagonalization of the matrix of coefficients. Instantons 
however do have a negative mode in their perturbation spectrum, as they describe unstable states. 
Since at any point relaxation amounts to solving the linearized equations for the correction around 
the present value of the solution, a relaxation procedure is bound not to converge. To take care of 
this problem what we did was to fix the value of the field at the center of the instanton, ip,,. (What 
this means in practice is that after each relaxation step we multiply the whole solution by a constant 
so that the center remains at the desired value without breaking the continuity of the relaxation 
process around it). The problem is then equivalent to that of a membrane with the center and the 
boundaries fixed to a certain value and subject to a potential, which should make it stable unless 
we have choosen too high a value at the center. In that case, the whole membrane will want to 
jump over the potential barrier (i.e., in the language of Ref. 7 we will be overshooting). Thus, 
if the value for (pc is too high the solution will still diverge, but if it is below critical relaxation 
should converge. At the boundary between these two cases should lie then the true value of the 
instanton at its center, ipTci ^^'^ fixing to this value should make the relaxation of the instanton 
a convergent process. 

We tried this procedure with the zero and high temperature solutions, where we could get 
independent figures for ipTc-i with a much higher precision from the shooting routine (since in both 
limits the equation of motion turns into a simpler differential equation). In practice we found that 
if the value chosen for was too high by an amount say A (so that ipc — ipTc'^^) the solution 
would start first by converging up to an overall accuracy level roughly of the order of A, and only 
after reaching that accuracy level would it start to diverge. Thus it usually didn't take long to 
produce a value for (p^ within a 0.1% difference from the ipxc value given by the shooting routine. 
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and a correspondingly low level of error for the whole solution. 



5 Results 

Using this algorithm then we aim at finding all the relevant finite temperature solutions of the 
equations of motion 

to elucidate whether the transition between the thermal hopping and the low temperature regimens 
is first or second order. We will do this for three different values of the parameter /. The first value 
that we will take for / will be chosen so as to check our procedure against Garriga's predictions for 
the thin wall case. Although in order to do this it would be desirable in principle to pick a value for 
/ as small as possible so that we get a wall as thin as possible, practical considerations limit how 
small a value for / we can choose. Since in the thin wall limit the wall radius is proportional to f~^, 
the smaller the value of / the larger the mesh that we would need to cover the whole bubble. But 
at the same time we need to have the mesh size small enough so that it can adequately represent 
the structure of the wall in a smooth way. Clearly then, if we do not want to have so many points 
in the mesh that the whole algorithm becomes unmanageable, there is a limit as to how small we 
can set the value of /. It turns out however that for our pourposes / = 0.25 is more than adequate. 

For such a value of /, we get a wall about half as thick as the interior of the bubble. Fig.5 
represents the shape of the potential for this case, an Fig.6 shows several finite temperature solu- 
tions. To compute the action of each solution we used a two dimensional version of the Romberg 
integration code that can be found in Ref. 14. When comparing the results obtained in this way 
with those provided by the shooting routine for zero temperature and the thermal hopping regimes, 
we consistently found a difference of about 0.2% of the total value. We feel that a conservative 
stimate for the error in the value of Se at each point should then be no higher than 0.5%. Fig.7 
gives us the Se vs. T dependence. From it we can deduce an approximate value for the temperature 
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Tc at which the transition takes place of Tc ~ 0.0476 in our adimensional units, which is in excellent 
agreement with the results from Garriga. Use of equation (24) above yields Tc ~ 0.0474, where 
we have used a value of R measured from the center of the bubble to the middle of the wall. The 
difference between the two values is well within our error stimates, and the transition is clearly seen 
to be first order. 

The reason behind such a good agreement with the thin wall results despite the fact that the 
wall is only half as thick as the bubble lies in Figs. 6b and 6c. In them we can see two complete 
periods of the solution for T — 0.91Tc and T — Tc respectively. We see that even at T = Tc the 
two spheres are far from touching each other, and due to the strong exponential decay outside the 
wall they almost do not interact with each other at all. In these conditions the ansatz used by 
Garriga for obtaining Eqn.23 (namely that the action of this solution is still the action for the zero 
temperature bounce) is obviously still valid, thus the agreement. Finally Fig.6d shows the thermal 
hopping solution that becomes dominant at temperatures T > Tc. 

The second value of / that we shall use will lay somewhere in the middle of the spectrum of 
possible values, at / = 0.55. Fig. 8 shows the shape of the potential for this case, clearly away from 
the region where the thin wall approximation is valid. Fig. 9a shows the zero temperature form of 
the solution (obviously already within the thick wall regime), and Figs. 9b and 9c show two complete 
periods for T — 0.9Tc and T — Tc. Here we can see some differences from the previous thin wall 
case: the two spheres collide and start merging into the wiggly cylinder structure, as shown more 
clearly in the aereal view depicted in Fig.9d, but at the critical temperature the solution (which in 
fact looks more like a "sierra" than like a cylinder) is still quite far away from the thermal hopping 
profile shown in Fig.9e. Finally in Fig. 10 we see that, although the transition is still first order it 
is obviously much less strongly so that in the previous case. The same considerations about errors 
apply in this case as above. 

The final value that we picked for / is very close to its upper possible hmit, / = 0.75 (for 
/ = 0.8 the double well structure of the potential totally disappears and we no longer have neither 
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a metastable state at ip- nor, consequently, any tunneling). The shape of the potential is again 
depicted in Fig. 11, where we see that in this case the field tunnels over a very shallow barrier. 
Fig. 12a shows the shape of the zero temperature bounce in this case, where we can see that the 
thick wall features of the solution have been accentuated. Fig. 12b shows again two complete 
periods at T = 0.6Tc. We see however in this case that the term proportinal to dip/dr has affected 
the sructure of the bounces, that no longer are spherically symmetric in this case, but stretched 
in the r direction. The peaks then start to interact "earlier" than in the previous cases, and 
at T = 0.6Tc they have already stablished the wiggly cylinder or "sierra" pattern that Linde 
advanced. This pattern continues to develop in Figs. 12c and d, where we can see how the solutions 
at T = 0.96Tc and T = 0.999Tc (the latter already inside the margin error for Tc itself) clearly and 
smoothly interpolate between the periodic array of bounces and the thermal hopping solution. The 
corresponding graph for the action is shown in Fig. 13. The figure shows a typically second order 
transition structure, with the curve for the action of the quantum regime merging smoothly up to 
the second derivative into the thermal activation curve. 

6 Conclusions 

In this paper we have addressed the issue of whether the transition from the thermal hopping regime 
to the quantum regime for bubble nucleation is first or second order, in the context of field theory 
</7^ phase transitions. Our main conclusion is that, whereas for most of the possible values of the 
parameter / that breaks the degeneracy of the two vacua the transition is first order, second order 
transitions will still occur when the potential barrier is sufficiently shallow. 

As mentioned in the introduction, any realistic model of the field phase transition will have to 
incorporate a temperature dependence in the potential, the specifics of which will depend on the 
particular values of the parameters chosen for the model. It seems clear from our results however 
that any model that does not settle sufficiently rapidly in the very shallow barrier regime will most 
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likely present a first order transtion between the thermal and the zero temperature regime, for only 
very shallow barrier penetration produces the very thick bubbles that seem to be needed in order to 
have smooth merging of the instantons and a second order transition. Indeed, were we to expect the 
results from the quantum mechanical case to carry over into the the field theory problem, one could 
have made the point that this actually had to be the case: since in the latter situation the point 
at which the field exits the barrier after tunneling is still quite far away from the stable vacuum 
and relatively close to the metastablc state, the field is not sensitive to the details of the potential 
beyond the barrier, and close to the true vacuum, while tunneling. In this case then one could 
approximate the potential by either a — c/?^ + </7^ or a — c/?^ + c/?^ expression. However, —(f + cf and 
—q^ + potentials are well known to produce second order transitions in quatum mechanics. 

The main practical difference between the two cases will lay in how abruptly the quantum 
tunneling rate will depart from the thermally induced one in the region close to the transition 
temperature, first order transitions marking a much sharper deviation from the thermal hopping 
regime than second order ones. In the context of condensed matter physics the same instanton 
techniques that we have applied here have yielded predictions for decay rates of metastable states 
via macroscopic quantum tunneling of the system at low temperatures, and about the nature of the 
transition from this to the thermal activation regime, that are amenable to experimental check (i.e. 
[15], [16-18] among others, and more recently [19]). In this context we are likely to observe very 
thick wall transitions, since in most of the situations of interest in experiments we either do not 
have a double well structure in the potential, or have rather shallow barriers. Indeed, second order 
transitions have been observed in Josephson junctions, and new experiments are being prepared 
to study the formation of bubbles of inversed magnetization in thin films under field inversion (a 
process completely analogous to the one that we have studied here only that in 2 + 1 dimensions) 
to check the predicted second order transition for that case. 
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8 Figure Captions 

Fig.l Motion of the particle in the inverted well: a tunneling trajectory between qi and q2 under the 
potential barrier V corresponds to a periodic motion in the inverted well —V with turning points 

Fig. 2 (a)Monotonic dependence of P on Euclidean Energy, (b)Euchdean Action vs. Temperature 
(~ (3~^) in this case 

Fig. 3 (a) Nonmonotonic dependence of /3 on Euclidean Energy, (b)Euchdean Action vs. Temper- 
ature (~ in this case 

Fig.4 Different shapes of the instanton as temperature increases: inside the shaded regions the 
field attains its true vacuum value, (a) Zero Temperature bounce, (b) Periodic array of bounces at 
T < 2R~^, (c) Wiggly cylinder, (d) Thermal hopping solution 

Fig. 5 Shape of the potential for f=0.25 

Fig. 6 Instantons describing the nucleation of the bubble for f=0.25. (a) T — 0, (h) T — 0.91 Tc, (c) 
T = r,, (d) T>T,. 

Fig. 7 Temperature dependence of the Euclidean Action: Sjnin{T)/S{0) versus T/Tf. for / = 0.25. 
Fig. 8 Shape of the potential for / = 0.55 

Fig. 9 Instantons describing the nucleation of the bubble for f=0.55. (a) T — 0, (b) T — 0.9 Tc, 
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(c) and (d) T = Te, (e) T>T,. 



Fig. 10 Temperature dependence of the Euclidean Action: Smin{T)/S{0) versus T/Tc for / = 0.55. 
Fig. 11 Shape of the potential for / = 0.75 

Fig. 12 Instantons describing the nucleation of the bubble for f=0.75. (a) T = 0, (b) T = 0.6 Tg, (c) 
T = 0.90 Te, (d) T = 0.96 Te, (e) T = 0.999 Te 

Fig. 13 Temperature dependence of the Euchdean Action: Sjnin{T)/S{0) versus T/Tc for / = 0.55. 
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